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Abstract 

This paper describes observations of a new phenomenon in evolved mass-losing 
AGB stars: the presence of two winds with different expansion velocities. CO(2-l) 
and CO (3-2) line emission was observed for 45 AGB stars at high velocity resolution 
and double winds found in 20% of the sample. Highly asymmetric lines were found 
in six other stars. The data tentatively suggest that double winds occur when the 
star undergoes a change (pulsational mode, chemical composition) and that the 
very narrow components represent the onset of a new phase of mass loss. 

1. Introduction 

Highly evolved stars on the asymptotic giant branch (AGB) and beyond are 
surrounded by extensive circumstellar envelopes produced by the copious mass loss 
which dominates the evolution at this stage. Observations of infrared continuum 
emission from circumstellar dust and millimeter wavelength line emission from cir- 
cumstellar molecules give the mass loss rates, wind outflow speeds, chemistry and 
structure of the winds and allow the mass loss history of the star to be traced. Sev- 
eral extensive data compilations which provide the empirical basis for investigating 
the mass loss phenomenon have been made (e.g. Loup et al. 1993), and recent 
reviews include those by Habing (1996) and Olofsson (1996, 1997a,b). 

Circumstellar molecules emit both maser and thermal radiation, and the latter 
is modeled with a fair degree of success by calculating the emission from a spheri- 
cally symmetric wind expanding at a uniform velocity at a uniform mass loss rate, 
with the effective envelope radius determined by molecular photodissociation by 
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the interstellar radiation field (Morris 1980; Mamon, Glassgold and Huggins 1988; 
Kastner 1992; Groenewegen 1994). This model predicts simple line shapes; optically 
thick emission lines are parabolic while optically thin lines are flat-topped, and the 
full width of the line at zero power is twice the terminal outflow speed of the wind. 

While this model (hereafter called the "steady wind" ) gives a decent first order 
representation of most circumstellar shells, recent, more detailed observations have 
found many shells which deviate from this simple description: 

1. Bipolar outflows. The shells around many mass- losing stars (e.g. VY CMa, 
AFGL 2688, and OH231.8+4.2) show flattened, bipolar structure on scales of ~ 10 16 
cm (e.g. Olofsson 1997a). 

2. Very fast molecular winds. The highest steady wind outflow speeds, up to 
~ 40 km s -1 , are found from a few supergiant stars such as IRC+10420 and VY 
CMa (Loup et al. 1993; Knapp et al. 1997b). A small number of stars have, as 
well as the steady wind component, molecular emission at much higher speeds, in 
excess of 200 km s -1 in the most extreme cases (Cernicharo et al. 1989; Gammie 
et al. 1989; Young et al. 1992; Jaminet et al. 1992; Alcolea et al. 1996; Knapp, 
Jorissen and Young 1997; Sahai and Nyman 1997). The fast winds do not have 
parabolic line shapes, and it is unlikely that they are produced by constant- velocity 
outflow. When observed with sufficient angular resolution, fast winds appear to be 
bipolar and to be flowing from the poles of dense disks or toroids (Neri et al. 1992; 
Yamamura et al. 1994). 

3. Interrupted mass loss. Olofsson et al. (1996) describe observations of detached 
molecular shells around several carbon stars, which appear to be produced on ther- 
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mal pulse timescales. 

4. Multiple- velocity and/or non-spherical outflows. Margulis et al. (1990) pointed 
out that not all molecular line profiles from evolved stars are adequately described 
by the steady wind model; some lines appear to be roughly gaussian or triangular 
in shape. Some few of these non-parabolic profiles have since been observed with 
sufficient angular and/or velocity resolution to demonstrate that they appear to 
have two or more components, centered at the stellar systemic velocity but with 
different outflow velocities. Among the stars showing this phenomenon are Mira 
(Planesas et al. 1990a,b), X Her (Kahane and Jura 1996) and RS Cnc (Jorissen 
and Knapp 1997). The profile shapes differ from those seen for stars with fast winds 
in that both components have roughly parabolic shapes and the outflow speeds lie 
within the range (3 - 30 km s -1 ) observed for the large majority of evolved stars. 
These stars thus appear to have two steady winds. 

Recent improvements in receiver sensitivity, frequency coverage, frequency res- 
olution and baseline stability make it feasible to carry out low noise, high velocity 
resolution observations of a significant sample of circumstellar envelopes to inves- 
tigate the frequency of occurrence of the multiple wind phenomenon and how it is 
related to the stellar chemistry, mass loss rate and evolutionary stage. This paper 
begins the discussion and analysis of a series of observations made to investigate 
multiple slow molecular winds in circumstellar envelopes. 

2. Observations 

This paper describes observations of the CO(2-l) and CO(3-2) lines made with 
the 10.4 m Robert B. Leighton telescope of the Caltech Submillimeter Observatory 
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on Mauna Kea, Hawai'i. The stars to be observed were selected from large surveys 
of CO emission (Margulis et al. 1990; Nyman et al. 1992; Loup et al. 1993; 
Knapp et al. 1997b) as having strong CO line emission, i.e. brighter than about 
1 K as observed with a 10 m telescope. We also confined the observations to stars 
north of about —30° because of the significant atmospheric opacity at submillimeter 
wavelengths. We can therefore expect to make observations with a high signal to 
noise ratio at high velocity resolution in a reasonable amount of observing time (2 
hours or less). The observations were mostly carried out in the period Dec 28 1996 
to Jan 2 1997, when we observed almost all CO-bright stars in the right ascension 
range available at that time of year, about 23 h to 10 h . The observed sample also 
contains several stars with weaker CO emission which were observed to fill empty 
places in the schedule and a small number of observations of stars in other parts 
of the sky obtained in the previous two years, during test and engineering time. 
The sample of stars described in this paper is thus very incomplete in terms of sky 
coverage, and we plan to complete these observations in future observing sessions. 

In all, high velocity resolution and (mostly) high signal-to-noise ratio observa- 
tions were obtained for 43 stars, which are listed in Table 1. Table 1 also contains 
data for two additional stars, both of which are known from previous observations 
to have double outflows: o Ceti (Planesas et al. 1990a,b; Knapp et al. 1997a), and 
IT 1 Gru (Sahai 1992). The CO(2-l) observations of RS Cnc are from Jorissen and 
Knapp (1997). 

Table 1 lists the star's most common name, the observed position, the spectral 
type, the variable type, the basic chemistry (O = oxygen rich, C = carbon star, S 
= S star), the period where known and a simple classification for the shape of the 
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CO line profile (see below). 

Most of the spectral types and chemistry classes are from Loup et al. (1993) 
or from the SIMBAD listings. Most of the variable types and periods are from the 
Catalogue of Variable Stars by Kholopov et al. (1985); the period of 55 d listed in 
that catalogue for EP Aqr is uncertain. The periods for several stars with thick 
circumstellar shells (R Scl, R Lep, AFGL 865, AFGL 971, IW Hya, IRC+10216 and 
IRC+40540) are given by Le Bertre (1992, 1993) and Cohen and Hitchon (1996). 
The period of the post-AGB carbon star HD 56126 is from Lebre et al. (1996). 
Several of the variable types in column 5 are from Kerschbaum and Hron (1992, 
1994); these authors show that the semi-regular SRa class (Kholopov et al. 1985) 
contains a mixture of Mira and SRb (hereafter "true" semiregular variables) and 
have reclassified several stars; R Scl from SRa to SRb and BK Vir and IT 1 Gru from 
Lb (cool irregular variable, see Kerschbaum et al. 1996) to SRb. 

The molecular line observations were made in the CO(2-l) or CO(3-2) line 
depending on the weather, which was variable during the observing run. The ob- 
servations were made using liquid-helium cooled SIS junction receivers with double- 
sideband system temperatures of ~ 100 K and 130 K at 230 and 345 GHz. The 
telescope half-power beamwidths are 30" and 20" at these frequencies, and the 
main-beam efficiencies are 76% and 65%. The receivers for both frequency bands 
are double sideband, and spectral lines in the image sideband, 3 GHz below the 
observed frequency for both receivers, were also detected. 

The spectral lines were observed by three acousto-optic spectrographs (AOS), 
with bandwidths of 1.5 GHz over 2048 channels, 500 MHz over 1024 channels and 50 
MHz at 1024 channels. The velocity resolution of the 50 MHz AOS is ~ 0.12 km s" 1 
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at 345 GHz and ~ 0.2 km s -1 at 230 GHz. The spectrometer frequency was 
calibrated using an internally generated frequency comb, and the velocity scale is 
corrected to the Local Standard of Rest (LSR). 

The observations were made by chopping between the star position and an 
adjacent sky position with the secondary mirror, using a chop throw of 120" in 
azimuth at a rate of 1 Hz. Pairs of chopped observations were made with the source 
placed alternately in each beam. The spectral baselines resulting from this proce- 
dure are linear to within the r.m.s. noise for all bandwidths with a small number 
of exceptions for which the 500 MHz data show small (0.005 - 0.01 K) baseline 
non-linearities. The CO emission from all of the observed stars is strong enough 
that it was used to measure the telescope pointing offsets before each observation. 

The temperature scale and the atmospheric opacity were measured by com- 
parison with a hot (room temperature) load. The line temperature was corrected 
for the main beam efficiency and the resulting scale is the Rayleigh- Jeans equivalent 
main-beam brightness temperature, Tmb, i-e. that measured by a perfect 10.4 m 
antenna above the atmosphere. The high velocity resolution CO line profiles are 
shown in Figure 1. 

Most of the line profiles are closely approximated by parabolic or flattened 
parabolic shapes (e.g. AFGL 865, W Aql). These line profiles are classified in 
Table 1 as parabolic, P. The CO(2-l) line profile for IRC+10216 shows the small 
blue-shifted feature which can be attributed to radiative transfer effects in a shell 
with a monotonically decreasing radial temperature gradient and a small amount of 
turbulence (Huggins and Healy 1986). Two of the profiles, those for AFGL 971 and 
V Cyg, are affected by galactic emission and hence their profile shapes are uncertain. 
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There is strong galactic emission near km s _1 in the 'on' and 'off' positions seen 
in the broad-band profile of IRC+60144 which may affect the observed line shape. 
All three line profiles are classified as parabolic. The observations of several other 
stars (R And, S CMi, S Vir, V CrB and R Ser) were made with insufficient signal to 
noise ratio to measure the profile shape, and are classified as parabolic by default. 
Uncertain classifications are enclosed in parentheses. 

Several line profiles are distinctly non-parabolic. Most striking are those show- 
ing two winds, i.e. two components of different width centered at the same veloc- 
ity; examples are IRC+50049, X Her and EP Aqr. These are classified as 'D', or 
double-wind profiles. Several other stars, e.g. R Scl and R Leo, have asymmetric 
line profiles, with one side of the profile being brighter than the other. These are 
classified as A (for asymmetric) Blue or Red, depending on whether the red or 
blue side of the profile is stronger. The possible origin of these line profile shapes 
will be discussed elsewhere. The stars with asymmetric profiles will be included in 
the group with parabolic line profiles in the discussion below, where stars in this 
combined group are compared with those with double winds. 

In almost all cases of stars with double winds, the narrow component has a 
similar width to the emission lines from interstellar molecular clouds - compare the 
observations of the contaminated profile of AFGL 871 and of the double- wind profile 
for EP Aqr in Figure 1. Are these double wind line profiles, then, simply an artifact 
of superimposed emission from molecular clouds? We believe not, for the following 
reasons. First, mapping observations of several of these stars: o Cet (Planesas et 
al. 1990a, b); RS Cnc (Neri et al. 1997); X Her (Kahane and Jura 1996); n 1 
Gru (Sahai 1992); and EP Aqr (Knapp et al. 1998) show that both the broad 
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and narrow components are co-located, centered on the star and of small angular 
extent. Second, the central velocities of both the broad and the narrow components 
are closely the same (see Table 2 below). We therefore have high confidence that 
both the broad and narrow components observed in the emission profiles represent 
circumstellar gas produced by mass loss. 

How well can the presence of double winds be found in these observations? 
A series of simulated CO(2-l) line profiles, with two components centered on the 
same velocity but with different outflow velocities, was constructed to investigate 
this question (for details see Section 4c below). In all cases, the fast wind component 
has an outflow speed Vf = 10 km s _1 and a mass loss rate of 3 x 10 -7 M© yr _1 . 
The model envelope was assumed to be at a distance of 200 pc and to be observed 
with the CSO 10.4 meter telescope; the resulting spectral line is parabolic with 
a peak temperature of Tf = 0.65 K. The models were calculated with a velocity 
resolution of 0.1 km s -1 , an r.m.s. noise of 0.4 K, and flat baselines. 

A second component was then added, with varying peak temperature T$ and 
outflow speed V s . These simulations showed that: 

(1) If V g /Vf > 0.7, the resulting line profile is sufficiently parabolic in shape that 
the presence of two winds cannot be distinguished. 

(2) If 0.7 < V s /Vf < 0.5, the presence of the two winds can be seen in a high- 
quality line profile provided T s ~ Tf. 

(3) If V s /Vf < 0.5, the slow wind can always be distinguished provided that T s is 
greater than five times the r.m.s. noise. 
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A selection of the simulated profiles is shown in Figure 2. 

The bottom line is that we can only reliably detect a double wind when there 
is a reasonable contrast between the fast and slow winds, i.e. when the outflow 
velocities differ by a factor of two or more. We therefore cannot use the data in 
this paper to rule out the possibility that all stars have multiple winds with similar 
expansion velocities. This caution applies to all of the discussion and conclusions 
in the remainder of this paper. It is of interest in this context that high spatial 
resolution far-infared maps of two of the carbon stars in Table 1, U Hya and Y 
CVn, find double/detached shells (Waters et al. 1994; Izumiura et al. 1996) but 
the global CO line profiles for these stars (Figure 1) show no discernable velocity 
structure. 

3. Line Profiles and Stellar Properties 
a. Chemistry 

Table 1 lists 45 stars: 26 are oxygen-rich, six are S stars and 13 are carbon 
stars. Six of the stars have asymmetric profiles, and all of these save R Scl are 
oxygen-rich. Nine of the stars have double profiles: three of these are S stars, 
one is a carbon star and five are oxygen stars; however, the presence of a narrow 
component for the carbon star IRC+60144 is somewhat uncertain (see Figure 1). 
The remainder have parabolic line profiles (and a few are observed with insufficient 
sensitivity to define the line shape well.) Thus, at least 38% of the oxygen stars, 
50% of the S stars and 15% of the carbon stars have non-parabolic line profiles. 
These results show the trend noted by Margulis et al. (1990); the line profiles from 
oxygen-rich circumstellar envelopes have more structure than do those for carbon- 
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rich envelopes. 

b. Pulsational Properties 

Figure 3 shows the histogram of the periods of stars with double and with 
parabolic line profiles; the former tend to have shorter periods. This difference is 
also apparent when variable type is considered; most (20) of the 26 Mira variables 
have parabolic profiles, while only half of the 15 SR variables do. Further, while 
the period distribution of the whole sample of stars has a peak at around 300 days, 
this peak is far more pronounced in the stars with double winds. 

c. Infrared Colors 

Figure 4 shows the [12] - [25], [25] - [60] and K - [60] infrared colors for the stars 
in Table 1, showing the dependence both on profile shape and on chemistry. Here, 
the colors are referred to those of a m star, using the zero point flux densities in the 
IRAS Explanatory Supplement (1988). The 12|Um, 25/um and 60/um flux densities 
from which these colors are calculated are from the IRAS Point Source Catalogue 
(1988) except for the S stars, whose flux densities are given by Jorissen and Knapp 
(1997). The K band magnitudes are taken from the compilation by Gezari et al. 
(1993), from Le Bertre (1992, 1993) and from Kerschbaum and Hron (1994). There 
may be large uncertainties in these colors, especially in K - [60], because the flux 
densities are not observed at the same point in the star's variability cycle (typically, 
AK ~ l m , Le Bertre 1992). Figures 4a and 4c show [25] - [60] vs [12] - [25] (cf. van 
der Veen and Habing 1988) and K - [60] versus [12] - [25], respectively, showing the 
colors as a function of profile shape, while Figures 4b and 4d show these relationships 
as a function of envelope chemistry. Generally, the higher the mass loss rate the 
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redder the infrared colors (van der Veen and Habing 1988). Figure 4 shows that 
stars with double winds are not present in the regions of these color-color diagrams 
occupied by stars with high mass loss rates; that this is not simply an artifact of 
the dependence of infrared color on circumstellar chemistry is shown by comparison 
with Figures 4b and 4d. 

d. Summary 

Nine of the 45 evolved stars in the sample in Table 1 have double winds, and 
only one of these is a carbon star. The stars with double winds have infrared colors 
which suggest that they have mass loss rates lying at the low end of those in the 
sample, and also have shorter periods than average. Further, three of the six S 
stars have double winds. There are two tentative conclusions from these properties: 
(1) it is unlikely that the double wind phenomenon is the result of binarity, since 
the frequency of double winds appears to depend on stellar chemistry and variable 
type, both of which are intrinsic to the star: and (2) they suggest that the double- 
wind phenomenon may be associated with period irregularities and occur when the 
chemistry and perhaps stellar pulsation and mass loss are undergoing a transition 
from one phase to another. 

4. Physical Properties of the Circumstellar Envelopes 
a. Kinematics 

The parameters of the spectral lines, measured from the high velocity res- 
olution profiles, are given in Table 2, which lists for each star the CO transition 
observed, the r.m.s. noise, the integrated CO line intensity Ico in K x km s _1 , 
and the peak line temperature T pea k, the central velocity V c and the terminal wind 
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outflow velocity V . The last three quantities were calculated for most of the line 
profiles by fitting a flattened parabolic line profile of the form 

T(V) =T peak (l-exp[-a(l-(^-^) 2 )])/(l-e- Q ) (1) 

» o 

(Knapp and Morris 1985) to the data. This was also done for the asymmetric 
line profiles such as that for R Leo to yield average quantities. The line profiles 
for IRC+60144, AFGL 971 and V Cyg are partly contaminated by galactic emis- 
sion. Equation 1 was fit to the uncontaminated regions of the profile, and the inte- 
grated line intensity Ico was estimated assuming that the line shape is parabolic: 
Ico = 4/3T peak V . These values are enclosed in parentheses in Table 2. 

The line profiles for stars with two components were modeled by a superpo- 
sition of two parabolas. No fit was made for the CO (2-1) profile for R Hya; this 
profile is very asymmetric. CO(2-l) data for RS Cnc from Jorissen and Knapp 
(1997) and CO(3-2) data for o Cet from Knapp et al. (1997a) were also analyzed. 
The profile fits for the stars with double winds are shown in Figure 5. Table 2 
lists the parameters for the broad and narrow components separately. As Figure 5 
shows, these fits are adequate for most of the stars. The narrow components for 
several stars appear to be more gaussian than parabolic in shape. Two stars, o Cet 
and R Hya, have broad components whose redshifted emission is stronger than the 
blueshifted emission; the profiles for the remaining stars are roughly symmetric. 

The central velocities V c and outflow velocities V Q measured by the CO(2-l) 
and CO (3-2) lines agree well for the stars for which both lines have been observed. 
The central velocities of the broad and narrow components are the same to within 
the errors for all stars, with < V c (broad) — V c (narrow) > = —0.2 km s -1 . Thus 
both flows are centered on the star (and among other things this means that it 
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is very unlikely that the narrow components are due to a chance superposition of 
interstellar emission on the stellar line profile) . Figure 6 shows the histograms of the 
outflow velocities for stars with double winds: the outflow velocities of the broader 
component are shown in Figure 6a and of the narrower component in Figure 6b. 
For comparison, the histogram of outflow speeds for the stars in Table 2 with single 
winds is shown in Figure 6c. Figure 6c appears to contain data for a representative 
sample of stars; the distribution is very similar to that found for other samples (see, 
e.g. Zuckerman, Dyck and Claussen 1986). The distributions of outflow velocities 
in stars with double winds differ from those in stars with single winds (although 
this difference is subject to the difficulties in distinguishing two components when 
the outflow velocities have similar values - see the discussion in Section 2). The 
broad components have somewhat larger outflow velocities than typical, but much 
more striking is the difference between Figures 6b and 6c; the narrow components 
have much smaller outflow velocities than average. Figure 6 suggests that the 
really unusual property of most of the stars with double winds is the presence of a 
component with a very small outflow speed. (The exceptions are IRC+60144 with 
V s = 5.6 km s" 1 and V f = 20.5 km s _1 , and II 1 Gru, with V s = 12 km s" 1 
and Vf = 27 km s -1 .) Indeed, Figure 6 suggests that there are essentially no stars 
with single winds with outflow speeds less than about 3 km s -1 , and that these low 
velocity flows appear only in stars with double winds. 

Figure 6c also shows the distributions of the wind speeds in stars with double 
winds, i.e. the sum of the data from Figures 6a and 6b. The distributions for stars 
with single winds and with double winds are quite different. This implies that there 
is indeed an intrinsic difference between the envelopes of stars whose line profiles 
show one component and those which show two components. 
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Figure 7 shows the velocity of the slow wind versus that of the fast wind for 
the nine stars with double winds. With the exception of EP Aqr, the slow and fast 
wind velocities are proportional to each other. 

b. Distances 

In the next section, we model the circumstellar CO emission to derive values 
for the mass loss rates. A prerequisite for this modeling is a measurement or estimate 
of the stellar distance, and this situation has changed dramatically in recent months 
with the publication of the Hipparcos catalogue (Perryman et al. 1997 - see also 
van Leeuwen et al. 1997). Thirty-three of the stars in Table 1 are in the Hipparcos 
Main Catalogue, and 23 of these have parallaxes measured at greater than the 2a 
level. The distances to these stars are taken directly from the Hipparcos catalogue. 

Other methods for estimating distances to AGB stars include using a period 
- My or period - Mk relationship; assuming a constant absolute K magnitude; or 
assuming a constant bolometric luminosity (see, e.g. Claussen et al. 1987; Jura 
1988; Young 1995; Knapp et al. 1997b). These methods can be re-evaluated using 
the Hipparcos distances, and a preliminary discussion is given below. 

Several studies (e.g. Claussen et al. 1987, Jura 1988, Jorissen and Knapp 
1997) have assumed the value of Mk = —8.1 m found for carbon stars in the Mag- 
ellanic Clouds (Frogel et al. 1980). The values of M K for the stars in Table 1 with 
parallaxes, regardless of chemistry, are shown in Figure 8a. These were calculated 
from the K magnitudes in the compilation of Gezari et al. (1993) and parallaxes 
from the Hipparcos catalogue. Apart from the supergiant a Ori (M K = — 9.6 m ) the 
stars appear to have a mean absolute K magnitude about l m fainter than — 8.1 m . 
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However, all of these stars are losing mass and so have significant circumstellar 
extinction. The approximate K extinction was found from: 

A K = 3xl0 7 MtMoyr" 1 ) V'^km s" 1 ) + l(T 4 D(pc) (2) 

The first term is the circumstellar extinction and the second the interstellar extinc- 
tion (with D the distance in pc). These quantities were found by assuming a gas to 
extinction ratio Nh/Av = 2 x 10 21 cm -2 mag -1 for both the circumstellar and 
interstellar gas, and an extinction curve with R = Av/Eb-v = 3.1, which gives 
Ay/Ax = 8.9. The circumstellar extinction was calculated assuming constant 
mass loss rate and outflow velocity, so that n(r) ~ r -2 , and an inner radius for the 
dust shell of 5 x 10 13 cm. The extinction corrections for the stars in Figure 8a were 
calculated from the mass loss rates given by the present CO observations (see below) 
and the resulting distribution of Mk is shown in Figure 8b. It should be empha- 
sized that this calculation is very uncertain; the differences between circumstellar 
and interstellar matter, and between envelopes of different chemical composition, 
have not been considered; there is growing evidence (see, e.g., recent papers by 
Skinner et al. 1997 and Whitelock et al. 1997) that AGB mass loss is anisotropic; 
and the apparent K magnitudes are uncertain because of variability (e.g. Le Bertre 
1992). It is clear from Figure 8b that the extinction has been much overestimated 
for several stars, but the results do support the use of a single absolute K magni- 
tude, consistent with the Magellanic Cloud carbon-star value of — 8.1 m , for all stars, 
and show that circumstellar extinction is often significant. 

For stars in the present study which do not have Hipparcos parallaxes, the 
distances were found using one or more of the above methods. Stars whose distances 
were found using Mk = — 8.1 m were modeled in an iterative process; a first distance 
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estimate was made, the mass loss rate was calculated, the circumstellar extinction 
was calculated, and the distance re-estimated. In all cases, the distance estimates 
were checked to be sure that they are not discrepant with the lower limits given by 
Hipparcos. These distance estimates are preliminary, and should be refined using a 
larger sample of stars; the absolute magnitudes of carbon and oxygen stars should 
also be further investigated. Distances were also estimated in some cases using 
Lboi = 10 4 Lq. The only serious inconsistency was found for R Hya, whose "K 
magnitude" distance (115 pc) is significantly smaller that the limit of > 200-300 
pc from the Hipparcos parallax - see also van Leeuwen et al. 1997b). No distance 
estimate was made for the post-AGB star HD 56126. 

c. Mass Loss Rates 

The mass loss rates were calculated using the model of Morris (1980) and 
Knapp and Morris (1985) to calculate the emission from a model spherical envelope 
produced by constant mass loss at constant outflow speed. The envelope radius 
was assumed to be truncated by photodissociation of CO, using the rates given 
by Mamon, Glassgold and Huggins (1988). The CO/H 2 abundance was taken as 
5 x 10~ 4 for oxygen stars, 6.5 x 10 -4 for S stars and 10~ 3 for carbon stars (see 
Knapp et al. 1997b). The slow and fast components for stars with double winds 
were treated separately. The results are given in Table 3, which gives for each 
star the distance, the chemistry, the assumed outflow velocity V Q , the resulting 
mass loss rate M and the source of the distance. These mass loss rates generally 
agree reasonably well with previous estimates - that for IRC+10216, for example, 
is about a factor of 2 lower than the value found by Crosas and Menten (1997) after 
the different model abundances are taken into account. 
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Figure 9 shows histograms of the mass loss rates for the separate double-wind 
components: the broad component (Figure 9a) and the narrow component (Figure 
9b). For comparison, the distribution for the stars in Table 3 with single winds is 
shown in Figure 9c. The broad components have a distribution similar to that of 
the stars with single winds; the narrow components have mass loss rates which are 
much smaller than typical. 

Several authors (e.g. Netzer and Elitzur 1993; Young 1995) have demonstrated 
an empirical relationship between mass loss rate M and outflow speed V Q for winds 
from AGB stars which provides a useful method for displaying and comparing data 
sets. Such a relationship is expected for radiationally-driven dusty winds, since the 
grain-gas coupling becomes less effective for lower mass loss rates (Gilman 1972; 
Goldreich and Scoville 1976; Netzer and Elitzur 1993; Ivezic and Elitzur 1995; 
Crosas and Menten 1997). 

The mass loss rate and outflow speed are plotted against each other in Figures 
10 - 12; these figures show the same data subdivided in several different ways. 
Figure 10 shows the plot for stars with distances measured by Hipparcos and those 
measured by other means. Both samples show a strong correspondence between M 
and V , but Figure 10 raises the interesting possibility that the distances to the stars 
found to have the highest mass loss rates (most, though not all, are carbon stars) 
have been overestimated and hence the mass loss rates have been overestimated: 
these stars, with M > 3 x 10~ 6 M , do not have Hipparcos parallaxes (they are 
optically too faint) and lie above the relationship found for the stars which do. 
This should be further investigated with a larger sample. The stars with Hipparcos 
distances show M ~ V^, a shallower relationship than found by Young (1995). 
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Figure 11 shows the data broken down into single winds, the slow components 
of double winds and the fast components. All types of component appear to follow 
the same relationship between M and V . Figure 12 shows the data divided by 
stellar chemistry; all of the stars appear to follow the same relationship. 

5. Discussion 

The results presented in the previous sections of this paper show that the 
phenomenon of evolved stars with double winds is quite common: at least 20% 
of the stars observed with high velocity resolution have two components centered 
at the same velocity. The present data show that the phenomenon is particularly 
common in S stars but this remains to be confirmed by observations of a larger 
sample. In addition, several stars were found to have strongly asymmetric profiles. 
These phenomena are almost absent in carbon stars. The tendency of this sort 
of velocity structure to be found in semiregular variable stars, or stars with short 
periods, as well as in S stars, suggests that double winds accompany changes in the 
stellar properties (pulsation mode, period etc) (see also Olofsson 1997a; Jorissen 
and Knapp 1997). 

If all stars go through one or more stages in which winds at two different 
velocities are produced, the present data show that they spend at least 20% of their 
AGB evolution in such a phase. The total lifetime of stars in the mass-losing phase 
on the AGB has been estimated at ~ 2 x 10 5 years (e.g. Young et al. 1993a,b), 
so that stars would spend ~4x 10 4 years in a phase or phases in which two winds 
are produced. 

Studies of Mira variables have found a correlation between pulsation period 
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and mass loss rate (e.g. Jura 1988; van der Veen and Habing 1988), and between 
period and wind outflow speed (Dickinson and Chaisson 1973), demonstrating the 
role of pulsation in mass loss. Semi-regular variables, which generally have shorter 
pulsation periods, do not show a period-velocity correlation (Dickinson and Dinger 
1982). Figures 13 and 14 show these relationships for stars in the present sample 
whose periods are known (Table 1), sorted by variable type. The period - mass loss 
rate correlation for Mira variables in Figure 13 and the period - outflow velocity 
correlation for Mira variables in Figure 14 are fairly well defined, but are not followed 
by the semi-regular variables, as found by Dickinson and Dinger (1982). The data for 
Miras with double winds shown in these figures suggest that the considerable scatter 
in these relationships may be intrinsic, and that the stellar properties (luminosity, 
period etc.) do not uniquely determine the mass loss rate. 

As the data discussed in this paper have shown, it is the slow winds which 
are unusual; they have outflow speeds and mass loss rates far lower than typical 
for evolved stars. Further, as Figure 7, 13 and 14 show, both the slow and fast 
winds are related to each other and to the stellar properties. Most of the stars with 
double winds have fairly low mass loss rates. As a result, the CO-emitting extent 
of the circumstellar envelope, determined by photodissociation, is small. Table 4 
lists the photodissociation radii R p for both the slow and fast winds for the stars 
with double winds; these were calculated using the results of Mamon, Glassgold 
and Huggins (1988) as the radii at which the relative CO abundance has decreased 
to half of its original value. Table 4 also gives the crossing time T = R p /V for 
both fast and slow winds. If the double-line profiles are due to the cessation of 
mass loss at one outflow velocity/mass loss rate and its resumption at a different 
outflow velocity /mass loss rate, these times are upper limits to the time interval 
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which can elapse between the cessation of mass loss and its resumption. For several 
stars (e.g. RS Cnc) this interval is < 1000 years, far shorter than the approximately 
5000 years observed for some carbon stars with detached circumstellar envelopes 
(e.g. Olofsson et al. 1996; Lindqvist et al. 1996). 

However, the correspondence between the presence of discernable double winds 
and stellar period, variability, type and mass loss rate suggest that the phenomenon 
only occurs at particular stages in a star's evolution. By analogy with carbon stars 
with detached envelopes, where the observations show that mass loss has resumed 
at a much lower mass loss rate and outflow speed (Olofsson et al. 1996; Lindqvist 
et al. 1996), and because winds with very low velocities (V Q < 3 km s _1 ) seem 
to occur only in stars with double winds, we tentatively conclude that oxygen and 
S stars also undergo episodes of interrupted mass loss, and that the slow winds 
represent its resumption. 

The present observations add to the growing body of data which show that 
the kinematics of circumstellar envelopes are a lot more complex that they appear 
at first, that changes in the stellar variability modes may produce abrupt changes 
in the mass loss rate, and that these changes are fairly common. Future papers 
in this series will discuss further survey observations, the systematics of slow and 
fast winds, asymmetric profiles, the relationship to the detached shells seen in some 
carbon stars (Olofsson et al. 1996) and individual stars in detail. 

6. Conclusions 

The CO(2-l) or CO(3-2) line profiles from 43 nearby AGB stars with bright 
CO emission were measured with high velocity resolution and sensitivity, and com- 
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bined with observations of two more such stars from the literature. We found: 

1. At least a third of the line profiles are strongly non-parabolic in shape. Six stars 
have asymmetric lines with bright spikes at the red or blue extrema of the profile 
and nine have double profiles, i.e. two roughly parabolic lines of different width 
centered at the same velocity. 

2. The stars with non-parabolic profiles have lower pulsation periods, bluer infrared 
colors (indicating lower mass loss rates) and are much more likely to be semiregular 
variables, than most of the stars in the sample. Only two of the carbon stars (15%) 
have non-parabolic profiles, half of the S stars do, and 40% of the oxygen stars. 
These results suggest that complex winds, whose presence is shown by complex line 
profiles, occur when the star is undergoing a change in luminosity, pulsation mode 
or chemistry. 

3. The broad and narrow components of the stars with double winds, analyzed 
separately, follow the same mass- loss-rate/period/outflow velocity relationships as 
do the stars with parabolic profiles. The broad components have mass loss rates and 
outflow velocities typical of mass losing AGB stars, but the narrow components have 
far lower velocities and mass loss rates than are seen in any other sample of evolved 
mass-losing AGB stars. The narrow components may be due to the resumption of 
mass loss after it has been stopped by some change in the stellar properties. 

4. Comparison of the line widths of 29 SiO(8 — 7) (for oxygen stars) or CS(7-6) 
(for carbon stars) with the CO(3-2) line widths shows that the SiO/CS lines are 
systematically narrower than the CO line for stars with V Q < 10 km s _1 ; above 
this velocity the line widths are roughly equal. These data suggest that the winds in 
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stars with low mass loss rates do not reach their terminal velocities until relatively 
large distances, several xlO 15 cm. 

5. A search for very broad line wings in this sample of stars, like those associated 
with post AGB stars (for example AFGL 618), proved negative except (tentatively) 
in two stars, where it need further observations. The phenomenon of fast winds 
appears to be confined to stars which are evolving away from the AGB. 
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Appendix A. Search for Fast (V > 50 km s -1 ) Molecular Winds 

A small number of evolved stars have very fast (V Q > 50 km s _1 ) molecular 
winds in addition to the slow (V Q < 20 km s -1 ) 'steady' winds. These winds, some 
with speeds in excess of 200 km s _1 , have, to date, been found in stars which are 
evolving away from the AGB. There is evidence that the formation of fast winds 
occurs after that of slow winds, and that these fast winds are bipolar (Morris et al. 
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1987; Cernicharo et al. 1989; Gammie et al. 1989; Bachiller et al. 1991; Neri et al. 
1992; Young et al. 1992; Jaminet et al. 1992; Bujarrabal et al. 1994; Alcolea et al. 
1996; Knapp, Jorissen and Young 1997). 

The observations described in the present paper are sensitive enough, and 
almost always have sufficiently flat baselines, that the line profiles measured with 
the 500 MHz AOS provide useful information on the presence or absence of very 
fast winds. As a guide, CSO observations of AFGL 618 find that the ratio of the 
line brightness temperatures of the fast and 'steady' winds is —0.04 for the CO(2-l) 
line and 0.09 for the CO(3-2) line (Gammie et al. 1989), where the temperature of 
the fast wind emission is measured at ±50 km s _1 from the central velocity. (The 
difference in these ratios is probably due to the fast wind being confined to the inner 
regions of the CRL 618 circumstellar envelope while the envelope produced by the 
steady wind is more extended and is partly resolved). CSO observations of AFGL 
2688 give ratios of about -0.12 and 0.24 (Young et al. 1992), and of V Hya -0.07 
and 0.06 (Knapp, Jorissen and Young 1997). 

However, a word of caution is in order here - the observations in this paper give 
information only on the fast wind emission from evolved stars. Sahai and Nyman 
(1997) describe their beautiful observations of the "Boomerang" bipolar nebula in 
the CO (1-0) line where a fast wind is seen in absorption against the microwave 
background; the wind is presumed to be cooled by rapid adiabatic expansion to a 
temperature < 3 K. Gas at this temperature has negligible excitation of the CO 
rotational lines so is not detectable in absorption or emission in any line except 
CO (1-0). 

The 500 MHz filter bank observations of the stars listed in Table 1 were 
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reduced and calibrated as described in Section 2, and the reduced profiles were 
binned to a resolution of ~ 7 km s _1 and examined to see if high- velocity emission 
was present. The results are given in Table 5. The only two stars with any hint 
of such emission are R Lep and RX Boo, whose CO line profiles are shown in 
Figure 15; these stars should be re-observed, since this structure could be due to 
weak emission lines from other molecules in the circumstellar envelope. The CO (2- 
1) profile of IRC+10216 has weak high-velocity wings, with T(wing)/T(peak) ~ 
0.0048 - see Table 5. However, this feature is almost certainly instrumental; the 
frequency calibration comb, as measured by the AOS, shows wings at about this 
level (Table 5). This represents the real limit to the measurable brightness of broad- 
wing emission. 

Otherwise, no high-velocity emission is present for any of the other stars. 
Table 5 lists the 5cr upper limit on the brightness temperature of fast wind emission 
and its ratio to the peak temperature of the steady wind, also taken from the binned 
profiles. Broad-band data are not available for a small number of stars because the 
baselines are non-linear. The ratios in Table 5 are less than typically observed in 
stars with fast winds (see above) for the large majority of the stars. Almost all of 
these stars have late spectral types, showing that they are still on the AGB. The 
single exception is HD 56126, whose spectral type is F5Iab and which is probably a 
post- AGB star (van Winckel et al. 1996; Justtanont et al. 1996; van Winckel 1997). 
These data provide additional evidence (subject to the possible existence of cold, 
fast winds discussed above) that the stage when fast molecular winds are produced 
is short-lived, and occurs at an evolutionary stage beyond the AGB. 
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Appendix B. Comparison of SiO/CS and CO Outflow Velocities 

The broad band (500 MHz) CO (3-2) observations detected emission from the 
CS(7-6) line from carbon stars or the 29 SiO(8 — 7) line from oxygen stars in the 
image sideband: the rest frequencies of these lines are 342.883 GHz and 342.979 GHz 
respectively. Five of the carbon stars and seven of the oxygen stars have emission 
of sufficient strength in these lines to allow a comparison of the line widths. The 
outflow speeds V were calculated by fitting a parabolic line to the CS/SiO and 
CO line profiles. Note that at the approximately 1 km s -1 resolution of the 500 
MHz AOS, any structure in the line profiles is not resolved and the fitted outflow 
velocities are likely to find an average of the slow and fast wind velocities. 

The results are listed in Table 6, which gives the values for V Q for the CO(3-2) 
and 29 SiO or CS lines calculated from the 500 MHz line profiles, plus the profile 
classifications and the CO outflow velocities from the 50 MHz observations (see 
Section 4). 

Comparison of the CO velocities measured with the 500 MHz and the 50 MHz 
AOS spectrometers shows the effect of velocity resolution; the lines are broadened 
by about 1 km s -1 . However, the high resolution CO lines are broader than the 
CS/ 29 SiO lines, showing that the velocity differences in Columns 2 and 3 of Table 
6 are real. Three of the stars in Table 6, RS Cnc, R Hya and X Her, have double 
winds; the data tentatively suggest that both the slow and fast winds contribute to 
the 29 SiO emission for RS Cnc and X Her, and that the 29 SiO emission arises only 
from the slow wind in R Hya. This question should be re-examined with higher 
velocity resolution observations of the CS and SiO line profiles. 
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Figure 16 shows the CS/SiO velocity versus the CO velocity for this sample of 
twelve stars. Overall, the CO lines are slightly broader than the CS/SiO lines, prob- 
ably due to radiative transfer effects in the more extended CO envelope. Addition- 
ally, the stars fall into two groups; three carbon stars with V Q (CO) = V (CS) and 
V > 10 km s _1 ; and the remaining nine stars in which V (CS, SiO) < V o (C0) 
and V < 10 km s _1 . The second group contains both carbon and oxygen stars 
and stars with both parabolic and double line profiles. 

CO is abundant and has effective self-shielding against photodestruction; fur- 
ther, its rotational levels are relatively easily excited. 29 SiO and CS have much 
lower abundances and also need higher densities to produce significant rotational 
line emission; accordingly, the CO emission preferentially measures gas at much 
larger distances from the star than do lines of less abundant molecules. In these 
stars, then, the dense gas in the inner regions of the envelope, which gives rise to 
the CS and 29 SiO lines, may not yet be accelerated to the terminal velocity. The 
results in Figure 16 suggest that the winds from these circumstellar envelopes do 
not reach their terminal outflow velocity until distances of a few xlO 15 cm from the 
star. 
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—02 32 32.6 


A /TT 

M7 


SKb 


(J 


1 r n 

150 


A/Ked 


TD I 1 , . 

K rlya 


1 O 0£? CO AO 

13 zd 58.48 


oo ni o A r 
—ZA 01 24.5 


A /T£? r ~ 

M6.5e 


Mira 


r\ 
U 


000 

388 




S Vir 


i o on oo o 

13 3D 23.2 


r\r* r-r* 101 

—06 56 18.1 


a /tt~ 

M7e 


Mira 


r\ 

y) 


O TO 

378 


( P ) 


T T 7 TT 

W Hya 


13 46 12.2 


—28 07 06.5 


A /TO 

M8e 


SKa 


/ \ 

u 


OOO 

382 


T ) 

r 


D V I ) 

Ka rJoo 


14 21 56.7 


1 O F. F. F. /lor 

+25 55 48.5 


A /TT r. 

1V17.5 


SKb 


(J 




1 > 

r 


S CrB 


1 r in 01 CO 

15 19 21.53 


1 1 00 A a r 

+31 32 46.5 


a /tt~ 

M7e 


Mira 


(J 


o£?n 

360 


r 


V UrrJ 


1 r /IT A A f\Q 

15 4( 44. U8 


1 on a 00 t 
+39 4o zz. ( 


r^ft 0^ 
Oo, 2e 


Mira 




358 


( P ) 


K oer 


1 r /to OO O 

15 48 Zo.z 


i 1 c 17 no t 
+ 15 1 f U2. ^ 


a /T£? r. ~ 

Mo.5e 


Mira 


u 


O r.^7 

35 1 




X Her 


16 01 08.8 


+47 22 35.8 


M6 


SKb 





95 


D 


RU Her 


16 08 08.6 


+25 11 59 


M7e 


Mira 





485 


P 


W Aql 


19 12 41.68 


-07 08 08.4 


S3, 9 


Mira 


S 


490 


P 


X Cyg 


19 48 38.48 


+32 47 11.8 


S7, 2e 


Mira 


S 


407 


P 


VCyg 


20 39 41.32 


+47 57 45.4 


C5,3e 


Mira 


c 


421 


P 


EP Aqr 


21 43 56.46 


-02 26 40.9 


M8 


SKb 





55 


D 
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Table 1, continued 

Star a(1950) 5(1950) Sp. Var Chem P(days) Profile 

II 1 Gru 22 19 41.13 -46 12 02.4 S5, 7 SRb S D 
IRC + 40540 23 32 01.3 +43 16 27 C8,3.5 C 628 P 

R Cas 23 55 51.7 +51 06 36.4 M7e Mira O 431 A/Red 

Notes to Table 1: 

1. Spectral types, variable types, chemistry and periods mostly from Loup et al. 
(1993) and from SIMBAD listings (periods primarily from Kholopov et al. 1985). 
Periods for R Scl, R Lep, CRL 865, and CRL 971 from Le Bertre (1992). Period for 
IRC+50096 from Jones et al. (1990). Period for HD 56126 from Lebre et al. (1996). 
Period for IW Hya from Le Bertre (1993). Period for IRC+10216; mean of periods 
from Le Bertre (1992) and Cohen and Hitchon (1996). Period of IRC+40540 from 
Cohen and Hitchon (1996). 

2. Mira (o Cet): CO profile classification from data given by Planesas et al. 
(1993a,b) and Knapp et al. (1997). IT 1 Gru: CO profile classification from data 
given by Sahai (1992). 

3. CO line profiles for CRL 971, IRC+60144 and V Cyg affected by Galactic 
emission. 



4. CO line profile classification: 



37 

P: parabolic. Classifications in parentheses are for profiles with low signal-to- noise 
ratio, significant Galactic contamination or uncertain classification. 

D: double; two distinct components centered at the same velocity with different 
outflow velocities. 

A: asymmetric, with Red or Blue emission wings. 
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Table 2. CO Line Parameters from High Velocity Resolution Profiles 



Star 


line 


r.m.s. 


Ico 






V c 


V 


IRC + 40004 


2 - 


1 


0.064 


14.0 ±2.6 


0.51 ±0. 


.08 


-20.5 ± 1.8 


19.6 ± 2.4 


R And 


3 - 


2 


0.46 


12.8 ±3.6 


1.07 ±0. 


.63 


-15.3 ± 5.4 


8.8 ±4.5 


R Scl 


2 - 


1 


0.017 


27.8 ±0.5 


1.01 ± 0. 


.10 


-19.0 ±0.9 


16.5 ± 1.1 


IRC + 50049 


2 - 


1 


0.055 


8.1 ± 1.0 


0.50 ±0. 


10 


-2.5 ±0.5 


9.5 ±0.5 












0.52 ± 0. 


,08 


-1.8 ± 0.2 


2.6 ± 0.3 


o Cet 


3 - 


2 


0.23 


76.2 ±3.0 


5.8 ± 1, 


.2 


±46.0 ±1.0 


6.7 ± 1.0 












9.0 ± 1. 


,4 


±46.6 ±0.2 


2.4 ± 0.4 


IRC + 50096 


3 - 


2 


0.21 


32.8 ±6.3 


1.64 ±0. 


25 


-16.4 ±1.6 


14.9 ± 2.7 


IRC + 60144 


2 - 


1 


0.050 


(14.3) 


0.45 ±0. 


.08 


-48.0 ±2.0 


20.5 ±2.1 












0.16 ± 0. 


06 


-45.0 ± 1.5 


5.6 ± 2.4 


TX Cam 


3 - 


2 


0.17 


32.7 ±6.5 


1.27 ±0. 


.21 


±11.2 ±2.1 


19.0 ±3.2 


R Lep 


2 - 


1 


0.046 


13.2 ± 1.6 


0.44 ± 0. 


.06 


±12.3 ± 1.2 


17.4 ± 1.5 


a Ori 


2 - 


1 


0.039 


8.2 ± 1.1 


0.32 ±0. 


.06 


±3.4 ± 1.2 


14.2 ± 1.2 


U Ori 


3 - 


2 


0.10 


7.2 ± 1.5 


0.71 ±0. 


,18 


-38.1 ± 1.3 


7.5 ± 1.2 


AFGL 865 


2 - 


1 


0.041 


16.5 ± 1.4 


0.81 ± 0, 


,06 


±43.1 ±0.7 


15.3 ±0.9 


AFGL 971 


3 - 


2 


0.14 


(15.4) 


0.77 ± 0. 


,19 


±3.3 ± 2.9 


14.9 ± 5.1 


GX Mon 


2 - 


1 


0.052 


22.9 ±2.0 


0.75 ±0. 


06 


-9.2 ± 1.0 


19.5 ± 1.2 


R Gem 


2 - 


1 


0.038 


2.68 ±0.84 


0.05 ±0. 


02 


-59.2 ±0.7 


11.0 ±1.1 












0.30 ±0. 


.05 


-59.0 ±0.6 


4.8 ±0.9 


HD 56126 


2 - 


1 


0.048 


9.2 ±0.9 


0.64 ±0. 


.07 


±73.0 ±0.8 


10.7 ± 1.1 


S CMi 


3 - 


2 


0.095 


1.5 ±0.6 


0.31 ±0. 


.08 


±51.5 ± 1.4 


3.3 ± 1.4 
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Table 2, continued 



Star 


line 


r.m.s. 


Ico 


RS Cnc 


2 - 


1 


0.074 


12.5 ±1.2 




3 - 


2 


0.12 


18.0 ±1.9 


R LMi 


2 - 


1 


0.025 


2.72 ±0.38 


IW Hya 


3 - 


2 


0.122 


12.3 ±3.4 


R Leo 


3 - 


2 


0.076 


21.9 ±1.2 


IRC ± 10216 


2 - 


1 


0.28 


424.2 ± 8.3 


U Hya 


2 - 


1 


0.044 


10.0 ±0.6 


R Crt 


2 - 


1 


0.050 


11.1 ± 1.1 


BK Vir 


2 - 


1 


0.039 


8.4 ±0.5 


Y UMa 


2 - 


1 


0.030 


3.2 ± 0.3 


Y CVn 


3 - 


2 


0.12 


10.7 ±1.9 


RT Vir 


2 - 


1 


0.023 


4.9 ±0.4 


SW Vir 


2 - 


1 


0.022 


11.7 ±0.4 


R Hya 


2 - 


1 


0.022 


4.9 ±0.5 




3 - 


2 


0.10 


22.2 ±2.2 


S Vir 


2 - 


1 


0.025 


0.64 ±0.25 


W Hya 


3 - 


2 


0.17 


28.3 ±2.7 



Tpeak 


V c 


Vo 


188 ±0.0.14 


±7.5 ±0.5 


8.0 ±0.5 


1.00 ±0.13 


±7.2 ±0.2 


2.4 ±0.3 


1.10 ±0.2 


±7.5 ±0.3 


8.0 ±0.6 


1.66 ±0.2 


±7.2 ±0.3 


2.8 ±0.4 


0.24 ±0.06 


0.0 ±0.6 


7.8 ±0.5 


0.61 ±0.18 


±40.4 ± 2.5 


13.6 ±4.1 


2.0 ±0.2 


-0.5 ±0.5 


6.8 ±0.7 


19.1 ±0.6 


-25.5 ±0.3 


14.6 ±0.3 


0.98 ±0.10 


-30.9 ± 0.4 


6.7 ±0.5 


59 ± 07 


-1-11 2 + 07 


10 8 ± 7 


0.29 ±0.05 


±17.1 ±0.8 


5.8 ± 1.1 


0.40 ±0.05 


19.0 ±0.5 


5.4 ±0.6 


0.88 ±0.17 


21.1 ±0.9 


7.8 ± 1.3 


0.37 ±0.03 


±17.5 ±0.5 


8.9 ±0.6 


0.89 ±0.11 


-11.1 ±0.6 


7.8 ±0.7 


1.2 ±0.3 


-10.0 ±0.8 


11.0 ± 1.2 


0.8 ±0.2 


-9.7 ±0.6 


5.0 ±0.7 


0.09 ±0.03 


±9.8 ±2.3 


5.2 ±2.2 


2.0 ±0.3 


±40.8 ± 0.6 


8.1 ±0.8 



40 



Table 2, continued 



Star 


line 


r.m.s. 


Ico 


Tpeak 


V c 


Vo 


RX Boo 


3- 


2 


0.15 


23.9 ±2.4 


1.65 ±0.21 


±1.5 ±0.7 


9.6 ± 1.0 


S CrB 


2 - 


1 


0.034 


2.53 ±0.51 


0.25 ±0.04 


±1.1 ±0.9 


7.4 ± 1.2 


V CrB 


3- 


2 


0.23 


4.3 ± 1.0 


0.39 ± 0.08 


-99.4 ±1.3 


6.3 ± 1.2 


R Ser 


3 - 


2 


0.17 


3.5 ±1.9 


0.52 ±0.04 


±31.6 ±0.5 


5.3 ±0.9 


X Her 


2 - 


1 


0.026 


6.4 ±0.4 


0.38 ± 0.05 


-72.8 ±0.8 


8.5 ±1.0 












0.45 ±0.07 


-73.2 ±0.4 


3.2 ±0.5 




3 - 


2 


0.071 


13.3 ±1.3 


0.7 ±0.1 


-73.2 ±0.5 


9.0 ± 1.0 












0.99 ±0.12 


-73.1 ±0.3 


3.5 ±1.4 


RU Her 


2 - 


1 


0.011 


2.3 ±0.2 


0.18 ±0.03 


-12.1 ±0.8 


9.4 ± 1.5 


W Aql 


2 - 


1 


0.030 


43.9 ±1.2 


1.63 ±0.09 


-24.1 ±0.6 


17.6 ±0.8 


X Cyg 


2 - 


1 


0.038 


28.8 ±0.7 


1.98 ±0.14 


±9.9 ±0.4 


8.9 ±0.5 


VCyg 


2 - 


1 


0.067 


(25.0) 


1.61 ±0.09 


±14.0 ±0.5 


11.8 ±0.6 


EP Aqr 


9 


i 
± 


0.035 


12.8 ±0.8 


0.70 ±0.10 


-33.0 ±0.5 


10.8 ±0.8 












1.36 ±0.21 


-33.8 ±0.2 


1.4 ±0.2 




3 - 


2 


0.10 


24.3 ±2.3 


1.30 ±0.3 


-33.0 ±0.5 


10.8 ± 1.0 












3.05 ±0.60 


-33.7 ±0.2 


1.5 ±0.3 


n 1 Gru 


2 - 


1 








-10 


27 














-13 


12 


tC + 40540 


3 - 


2 


0.17 


49.8 ±4.8 


2.61 ±0.26 


-16.8 ±0.9 


14.3 ± 1.4 


R Cas 


3 - 


2 


0.20 


51.8 ±4.8 


2.94 ± 0.34 


+24.9 ±0.9 


12.1 ± 1.2 
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Notes: 

1. IRC+60144, AFGL 971 and V Cyg: contaminated by galactic 

2. R Hya CO(2-l) - no fit possible 

3. n 1 Gru: data from Sahai (1992) observed with 15 m SEST 
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Table 3. Mass Loss Rates 



Star 


D(pc) 


Chem 


V G (km s" 1 ) 


M(M yr" 1 ) 


Dist. Method 


IRC + 40004 


850 


O 


19.6 


8.4 x 10" 6 


Mk, L bo i 


R And 


490 


S 


8.8 


9.4 x 10" 7 


Mk 


R Scl 


475 


c 


16.5 


2.7 x 10" 6 


Mk, Hipparcos 


IRC + 50049 


142 





9.5 


1.5 x 10" 7 


Hipparcos 








2.6 


2.7 x 10" 8 




o Cet 


128 





6.7 


4.4 x 10" 7 


Hipparcos 








2.4 


9.4 x 10" 8 




IRC + 50096 


650 


c 


14.9 


4.4 x 10" 6 


Mk, L b oi 


IRC + 60144 


850 


c 


20.5 


4.8 x 10" 6 


Mk, L b oi 








5.6 


2.4 x 10" 7 




TX Cam 


280 





19.0 


2.3 x 10" 6 


Mk, L b oi 


R Lep 


250 


c 


17.4 


5.2 x 10" 7 


Hipparcos 


a Ori 


131 





14.2 


3.1 x 10" 7 


Hipparcos 


U Ori 


300 





7.5 


2.9 x 10" 7 


M K , P — M K , Hipparcos 


AFGL 865 


1800 


c 


15.3 


2.2 x 10" 5 


Mk, Lboi 


AFGL 971 


1300 


c 


14.9 


7.6 x 10" 6 


Mk, L b oi 


GX Mon 


500 





19.5 


4.8 x 10" 6 


Mk, Lboi 


R Gem 


850 


s 


11.0 


4.1 x 10" 7 


M K 








4.8 


4.5 x 10" 7 




S CMi 


470 





3.3 


6.6 x 10" 8 


P-M k ;M k 
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Table 3, continued 



Star 


D(pc) 


Chem 


V (km s" 1 ) 


M(M yr" 1 ) 


Dist. Method 


RS Cnc 


122 


S 


8.0 


1.0 x 10" 7 


Hipparcos 








2.6 


2.3 x 10" 8 




R LMi 


270 





7.8 


1.6 x 10" 7 


P — Mk; Mk; Hipparcos 


IW Hya 


900 





13.6 


4.2 x 10" 6 


Mk; Lboi 


R Leo 


101 





6.8 


9.4 x 10" 8 


Hipparcos 


IRC + 10216 


150 


c 


14.6 


9.2 x 10" 6 


Mk, Lboi 


U Hya 


162 


c 


6.7 


9.2 x 10" 8 


Hipparcos 


R Crt 


300 





10.8 


6.9 x 10" 7 


M K , P — Mk, Hipparcos 


BK Vir 


176 





5.8 


6.9 x 10" 8 


Hipparcos 


Y UMa 


313 





5.4 


1.7 x 10" 7 


Hipparcos 


Y CVn 


218 


c 


7.8 


1.1 x 10" 7 


Hipparcos 


RT Vir 


138 





8.9 


1.1 x 10" 7 


Hipparcos 


SW Vir 


143 





7.8 


1.7 x 10" 7 


Hipparcos 


R Hya 


200 





11.0 


4.7 x 10" 7 


M K ; Hipparcos 








5.0 


7.3 x 10" 8 




S Vir 


420 





5.2 


8.4 x 10" 8 


Hipparcos 


W Hya 


115 





8.1 


1.7 x 10" 7 


Hipparcos 


RX Boo 


156 





9.6 


3.0 x 10" 7 


Hipparcos 


S CrB 


400 





7.4 


2.5 x 10" 7 


M K ; P — M K ; Hipparcos 


V CrB 


600 


c 


6.3 


2.5 x 10" 7 


Mk ; Hipparcos 
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Table 3, continued 



Star 


D(pc) 


Chem 


V (km s 


R Ser 


279 





5.3 


X Her 


138 





9.0 








3.4 


RU Her 


400 





9.4 


W Aql 


610 


s 


17.6 


X Cyg 


106 


s 


9.0 


V Cyg 


271 


c 


11.8 


EP Aqr 


135 





10.8 








1.4 


n 1 Gru 


153 


s 


27 








12 


IC + 40540 


700 


c 


14.3 


R Cas 


107 





12.1 



M(M© yr" 1 ) Dist. Method 



9.2 x 


10" 8 


Hipparcos 


1.1 X 


lO" 7 


Hipparcos 


3.4 x 


10" 8 




3.2 x 


lO" 7 


M K ; P — M K ; Hipparcos 


9.4 x 


10" 6 


M K 


2.3 x 


lO" 7 


Hipparcos 


9.4 x 


lO" 7 


Hipparcos 


2.3 x 


lO" 7 


Hipparcos 


1.7 x 


10" 8 




5.7 x 


lO" 7 


Hipparcos 


4.6 x 


lO" 7 




7.4 x 


10" 6 


Mk; Lboi 


5.2 x 


lO" 7 


Hipparcos 
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Table 4. Photodissociation Radii of Stars with Double Winds 



Star V (km s 1 ) R p (cm) T(years) 



o Cet 


6.7 


3.8 x 10 16 


1800 




2.4 


2.4 x 10 16 


3100 


IRC + 60144 


20.5 


1.4 x 10 17 


2170 




5.6 


1.7 x 10 17 


9700 


R Gem 


11.0 


3.5 x 10 16 


1000 




4.8 


5.2 x 10 16 


3390 


RS Cnc 


8.0 


1.8 x 10 16 


710 




2.6 


1.1 x 10 16 


1370 


R Hya 


11.0 


3.3 x 10 16 


950 




5.0 


6.0 x 10 16 


3750 


X Her 


9.0 


1.7 x 10 16 


580 




3.4 


1.1 x 10 16 


1050 


EP Aqr 


10.8 


2.3 x 10 16 


670 




1.4 


1.1 x 10 16 


2390 


II 1 Gru 


27 


3.6 x 10 16 


420 




12 


3.6 x 10 16 


950 
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Table 5. Search for Fast Molecular Winds in Observed Stars 



btar 


Line 


1 (rast J 


K 


IRC + 4UUU4 


2 — 1 


< U.U15 


< U.U28 


R A r,A 

rvAnQ 


q o 
o — z 


<. U.lz 


<. U.lz 


IRC + 5UU49 


z — 1 


< U.Ul 


< U.U2 


IRC + 5UU96 


o o 
o — z 


/AH/1 

< U.U4 


< U.U9 


T"D ^ i A A 

1KC + 6U144 


z — 1 


< U.Ulo 


< U.U3 


1a uam 


o o 
o — z 


< U.Uo 


^ r\ r\ A 

< U.U4 


K Lep 


z — 1 


(U.U2J 


(U.UOJ 


a On 


z — I 


< U.Ul 


< U.U3 


U On 


o o 
o — Z 


< U.U3 


< U.U4 


A.t OL 865 


z — 1 


< U.U15 


^ r\ no 

< U.Uz 


Ai* OL 971 


3 — 2 


< U.U4 


< U.Uo 


Ga Mon 


Z — 1 


< U.Ulo 


< U.Uz 


Iv Gem 


Z — 1 


< U.Ul 


< U.Uo 


rliJ OOlzO 


z — 1 


< U.Ul 


< U.Uz 


b CMl 


3 — 2 


< U.U3 


< U.2 


RS Cnc 


2-1 


< U.U3 


< U.U2 


RS Cnc 


3-2 


< U.U3 


< U.U2 


R LMi 


2 - 1 


< U.Ul 


< U.U4 


IW Hya 


3-2 


< U.U25 


< U.U4 


R Leo 


3-2 


< U.U25 


< U.Ul 


IRC + 10216 


2-1 


(o.i) 


(U.UU48) 
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Table 5, continued 



Star 

U Hya 
R Crt 
BK Vir 

Y UMa 

Y CVn 
RT Vir 
SW Vir 
R Hya 
R Hya 
S Vir 

W Hya 
RX Boo 
S CrB 

Y CrB 
R Ser 
X Her 
X Her 

RU Her 
W Aql 
X Cyg 
VCyg 



Line 

2-1 
2 - 1 
2 - 1 

2- 1 

3- 2 

2- 1 
2 - 1 
2 - 1 

3- 2 
2 - 1 
3-2 
3-2 
2 - 1 
3-2 
3-2 

2- 1 

3- 2 
2 - 1 
2-1 
2 - 1 
2-1 



T(fast) 

< 0.01 

< 0.01 

< 0.01 

< 0.015 

< 0.03 

< 0.02 

< 0.02 

< 0.02 

< 0.035 

< 0.01 

< 0.03 
(0.05) 

< 0.02 

< 0.06 

< 0.03 

< 0.03 

< 0.03 

< 0.01 

< 0.03 

< 0.02 

< 0.06 



R 

< 0.01 

< 0.02 

< 0.06 

< 0.06 

< 0.04 

< 0.062 

< 0.026 

< 0.06 

< 0.02 

< 0.13 

< 0.02 
(0.03) 

< 0.10 

< 0.14 

< 0.07 

< 0.06 

< 0.02 

< 0.06 

< 0.017 

< 0.009 

< 0.031 
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Table 5, continued 



Star Line T(fast) R 

EP Aqr 2-1 < 0.02 < 0.022 

IRC + 40540 3 - 2 < 0.03 < 0.01 

R Cas 3-2 < 0.04 < 0.01 

AOS comb 0.0061 



Note: (1) T(fast) = brightness temperature of fast wind. R = ratio of brightness 
temperatures of fast and steady wind. Doubtful values in parentheses. 

(2) Wings on the strong emission line from IRC+10216 are likely to be instrumental 
- compare the data for the AOS frequency calibration comb. 
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Table 6. Wind Velocities in the CO(3-2), CS(7-6) and 29 SiO(8 - 7) Lines 
a. Carbon Stars 

Star V o (C03-2) V G (CS7 - 6) Profile V o (C03-2) 

(50 MHz) 



IRC + 50096 


15.2 ±0.7 


13.7 ± 1.5 


P 


14.9 ±2.7 


AFGL971 


15.3 ±0.8 


14.4 ±4.4 


(P) 


14.9 ±5.1 


Y CVn 


8.7 ±0.6 


7.2 ±2.2 


P 


7.8 ± 1.3 


V CrB 


8.0 ± 1.5 


3.9 ±0.9 


(P) 


6.3 ± 1.2 


IRC ± 40540 


14.7 ±0.5 


15.7 ±1.6 


P 


14.3 ± 1.4 
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Table 6, continued 



b. Oxygen Stars 



Star V o (C03-2) V o ( 29 Si08 - 7) Profile V o (C03 - 2) 

(50 MHz) 



UOri 7.8 ±1.0 

RS Cnc 6.4 ± 1.7 



R Leo 
R Hya 

RX Boo 
X Her 

R Cas 



7.7 ±0.1 
8.9 ± 1.5 

10.0 ±0.5 
7.3 ± 1.9 

12.6 ±0.6 



5.6 ± 1.7 
4.2 ±0.2 

4.7 ±0.7 
3.7 ± 1.7 

7.1 ± 1.1 
5.7 ± 1.1 

5.2 ± 1.0 



P 
D 

A/red 
D 

P 
D 

A/red 



7.5 ±1.2 
8.0 ±0.6 
2.8 ±0.4 
6.8 ±0.7 

11.0 ± 1.2 
5.0 ±0.7 

9.6 ±1.0 
9.0 ±1.0 
3.5 ±1.4 

12.1 ± 1.2 
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Figure Captions 

Figure 1: CO line profiles observed at a velocity resolution of ~ 0.1 — 0.2 km s _1 
for 43 evolved stars. 

Figure 2: Simulated CO(2-l) line profiles as observed with the CSO 10.4 m telescope 
for stars with double winds. In all cases the faster wind has an outflow speed 
Vf = 10 km s _1 and a mass loss rate of 3 x 10~ 7 M yr _1 . The model star is at 
a distance of 200 pc. The simulated line profiles have a velocity resolution of 0.1 
km s _1 and an r.m.s. channel-to-channel noise of 0.04 K. The line profile of the 
fast wind is parabolic with a peak temperature of 0.65 K. (a) outflow speed of slow 
wind V s = 3 km s _1 , peak temperature T s = 0.2K (five times the r.m.s. noise), 
(b) V s = 6 km s _1 , T s = 0.6 K. (c) V s = 8 km s _1 , T s = 1.0 K. 

Figure 3: Period distribution for the evolved stars in Table 1. Shaded columns: 
stars with double molecular winds. 

Figure 4: [25] - [60] vs. [12] - [25] and K - [60] vs. [12] - [25] color-color diagrams 
for evolved stars. Figures 3a and 3c: colors as a function of line shape. Figures 3b 
and 3d: colors as a function of chemistry. 

Figure 5: CO profiles for stars with double winds showing profile fits (dashed lines). 

Figure 6: Histograms of wind outflow speeds for (a) the broad component of stars 
with double winds (b) the narrow components of stars with double winds and (c) 
stars with single winds. The summed distributions for the broad and narrow com- 
ponents is shown by dashed lines. 

Figure 7: Outflow speeds of slow versus fast winds for 9 stars with double winds. 
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Figure 8: Absolute K magnitudes for stars with Hipparcos parallaxes (a) observed 
(b) corrected for circumstellar and interstellar extinction. 

Figure 9: Distributions of mass loss rates for (a) broad components and (b) narrow 
components of stars with double winds: and (c) stars with single winds. 

Figure 10: Mass loss rate versus wind outflow speed. Filled symbols: stars with 
Hipparcos parallaxes. Open symbols: distances derived by other methods (see 
text). 

Figure 11: Mass loss rate versus wind outflow speed. Filled symbols: stars with 
double winds, broad components. Open symbols: stars with double winds, narrow 
components. Crosses: stars with single winds. 

Figure 12: Mass loss rate versus wind outflow speed. Filled symbols: carbon stars. 
Open symbols: oxygen stars. Crosses: S stars. 

Figure 13: Mass loss rate in M yr -1 versus pulsation period (days) for stars with 
measured periods (Table 1), sorted by variable type. Filled symbols: Mira variables. 
Open symbols: SRb variables. Crosses: SRa, SRc and unknown type. The data for 
a Ori (SRb: P = 2070 days) are not plotted. Vertical lines join the mass loss rates 
for stars with double winds. 

Figure 14: Wind outflow speed V Q in km s -1 versus pulsation period (days) for 
stars with measured periods (Table 1), sorted by variable type. Filled symbols: 
Mira variables. Open symbols: SRb variables. Crosses: SRa, SRc and unknown 
type. The data for a Ori (SRb: P = 2070 days) are not plotted. Vertical lines join 
the velocities for stars with double winds. 
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Figure 15: CO line profiles, binned to a resolution of ~ 7 km s -1 , for (a) R Lep and 
(b) RX Boo. 

Figure 16: Wind outflow speed measured from the CO (3-2) line versus that mea- 
sured from the 29 SiO(8 — 7) line (oxygen stars, open symbols) or the CS(7-6) line 
(carbon stars, filled symbols). 
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